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OVERVIEW 

We  have  developed  a  proposed  sensor  system  for  the  Topcat  anti-armor  system.  It 
analyzes  cross-sectional  measurements  of  the  underlying  terrain,  detecting  targets  by  their  dis¬ 
tinctive  shape.  Sparsely-sampled  point  from  a  single  scan  line  are  fed  as  input  to  a  Parzen  esti¬ 
mator,  which  makes  the  target/no-target  decision.  Morphological  operators  then  remove  spuri¬ 
ous  noise  points,  producing  large  blobs  over  the  detected  target  region.  Results  are  shown,  using 
the  GTVISIT  simulation  system  to  produce  test  data. 


RESULTS 

Topcat  is  a  proposed  smart  anti-armor  weapons  system.  It  would  be  a  missile  which  flew 
back  and  forth  over  the  battlefield  looking  for  targets  of  opportunity.  A  spinning  mm-wave  radar 
sensor  (or  LADAR  sensor)  in  its  nose  would  search  for  targets;  when  one  was  found,  the  missile 
would  fire  a  submunition  at  it.  We  have  developed  a  possible  algorithm  for  automatic  target 
recognition.  This  algorithm  has  been  tested  on  synthetic  data  produced  by  the  GTVISIT  simula¬ 
tion  system. 

Figure  1  shows  a  schematic  top  view  of  a  tank  and  a  tree.  A  scene  like  this  one  was  gen¬ 
erated  by  GTVISIT,  and  processed  into  simulated  Topcat  data  (this  involves  generating  a  new 
range  image  from  the  viewpoint  of  the  center  of  each  rotational  scan,  and  then  sampling  points 
lying  along  an  arc  in  this  image).  Figure  2  shows  the  range  profiles  associated  with  the  scans 
one  through  six  in  Figure  1. 

We  know  that  armored  vehicles  will  be  roughly  rectangular  in  cross-section.  Natural 
clutter  will  have  range  profiles  with  a  high  variance,  and  are  therefore  unlikely  to  exhibit  the 
characteristic  rectangular  shape  (or,  at  most,  only  in  very  isolated  cases).  There  are  three  general 
approaches  to  measuring  this  shape  and  making  the  associated  decisions: 

□  A  rule-based  system  that  analyzes  the  presence  or  absence  of  step  edges  and  straight  line 
segments. 

□  A  signal-processing  based  approach. 

□  A  statistical  pattern  recognition  type  of  classifier,  where  the  feature  vector  is  a  sequence  of 
range  (or  height)  measurements.  Since  the  amount  of  training  data  required  for  accurate 
classification  is  proportional  to  the  length  of  the  feature  vector,  this  should  be  kept  to  a 
minimum,  suggesting  sparse  sampling. 

We  elected  to  use  the  third  approach;  it  yielded  an  accurate  and  computationally  efficient 
method.  A  robust  non-parametric  classifier,  the  Parzen  classifier,  was  implemented.  Being  non- 
parametric,  this  classifier  does  not  entail  any  a  priori  assumptions  regarding  the  probabiliity  dis¬ 
tributions  of  the  height  data  generated  by  the  Topcat  sensor.  More  specifically,  unlike  the  more 
commonly  used  classifiers,  this  classifier  does  NOT  assume  that  that  the  height  values  for,  say. 
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the  trees  are  Gaussian. 

In  our  Parzen  classifier  for  single  scan  line  data,  we  use  N  adjacent  height  values  as 
features.  To  provide  further  explanation  regarding  classification,  we  will  denote  a  set  of  N  adja¬ 
cent  height  values  by  . ,/n ■  During  the  training  phase,  a  large  number  of  sets  of  N 

feature  values  are  collected  from  the  synthetic  data,  first  from  just  the  tanks  and  then  from  just 
the  background.  These  sets  of  N  values  correspond  to  single  points  in  an  N  dimensional  feature 
space.  During  the  training  phase,  let  the  vector  f-'2'1*  denote  i-th  such  point.  If  we  use,  say,  100() 
sets  of  these  N  values  derived  from  just  the  tanks,  then  the  following  function  constitutes  a  Par¬ 
zen  estimate  of  the  probability  density  function  corresponding  to  just  the  tanks: 

1000 

/W(0  =  A-xIMf-fT*) 

l 

where  p  is  the  probability  density  function  and  f  any  arbitrary  point  in  the  feature  space.  The 
function  h,  referred  to  as  a  kernel  function,  admits  many  different  possibilities,  such  as  rectangu¬ 
lar,  triangular,  Gaussian,  etc.,  each  different  choice  affecting  differently  the  rate  of  convergence 
of  the  estimated  density  function  to  its  true  value.  We  use  Gaussian  a  kernel  since  it  has  been 
demonstrated  to  yield  the  best  rates  of  convergence  for  a  given  number  of  samples  for  training. 
Note  that  just  because  we  are  using  a  Gaussian  kernel  for  h,  it  does  not  imply  that  the 
classification  technique  is  parametric;  the  shape  of  the  overall  density  function  can  still  be  arbi¬ 
trary.  The  factor  K  in  the  above  formula  is  a  normalization  constant;  its  purpose  is  to  make  cer¬ 
tain  that  the  area  under  the  probability  density  function  is  unity. 

We  can  similarly  collect,  say,  a  1000  sets  of  N  adjacent  heights  for  the  background  clutter 
such  as  trees.  If  we  represent  each  set  by  i‘[ees ,  we  can  form  the  following  Parzen  estimate  of  the 
probability  density  function  corresponding  to  the  background  trees: 

1000 

Ptreesit)  =  K*  £  h(f  -  t‘rees) 

1 

After  a  classifier  is  trained  by  estimating  ptank  and  p,reeS.  a  set  of  measurements  on  an  unknown 
target  can  then  be  classified  by  invoking  the  Bayes  rule.  In  its  simplest  implementation,  which 
assumes  that  the  prior  probabilities  associated  with  the  tanks  and  the  trees  are  identical,  the  set  of 
mesurements  would  correspond  to,  say,  a  tank  if  and  only  if  the  following  condition  is  satisfied: 

Ptank  (0  —  Ptrees(^) 

Note  that  this  approach  to  classification  captures  the  correlations,  at  multiple  orders,  between  the 
different  height  values.  It  seems  intuitively  plausible  that  the  adjacent  height  values  for  a  tank 
would  be  more  highly  correlated  than  the  adjacent  values  for  a  bush  or  a  tree 

The  algorithm  processes  the  data  as  follows.  First,  each  scan  line  is  median  filtered,  an 
optional  step  that  may  be  used  to  reduce  the  noise  present  in  a  real  sensor.  The  data  is  then  nor¬ 
malized,  to  make  the  decision  insensitive  to  missile  altitude.  Sparsely-sampled  feature  vectors 
are  then  extracted,  and  a  Parzen  classifier  is  applied.  Morphological  operators  then  reduce  the 
noise,  leaving  large  contiguous  regions  labeled  as  “target”  on  the  armored  vehicles.  Figure  3 
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shows  an  overview  of  our  method. 

To  address  some  of  these  steps  in  more  detail,  the  normalization  can  be  as  simple  as  sub¬ 
tracting  the  average  range  for  a  scan  line  from  all  points  on  it.  The  rectangular  shape  is  pro¬ 
duced  by  the  derivative,  not  the  absolute  range.  Therefore,  the  instantaneous  missile  altitude  has 
no  effect  on  the  decision  process  (except,  of  course,  for  any  effects  it  may  have  on  sensor  accu¬ 
racy  or  noise  characteristics). 

Once  the  data  has  been  median  filtered  and  normalized,  feature  vectors  must  be  extracted. 
Since  the  objects  of  interest  are  well  known,  the  overall  size  of  this  feature  vector  is  specified. 
For  instance,  a  Russian-built  T-72  measures  approximately  8.1  meters  comer-to-corner  across 
the  hull.  Therefore,  the  feature  vector  must  be  large  enough  that,  even  in  a  worst-case  scenario 
where  a  scan  line  hits  both  extreme  comers,  the  outermost  points  on  the  feature  vector  will  lie  on 
the  ground  and  the  interior  ones  on  the  vehicle  itself. 

With  the  available  information  on  the  sampling  rate,  spin  rate,  and  typical  altitude  of  the 
proposed  sensor  system,  it  appeared  that  the  following  would  be  a  useful  feature  vector.  Note 
that  r,  is  the  f  th  range  measurement  in  a  scan  line,  r  is  the  median-filtered  version  of  r,  and  X,  is 
the  feature  vector  for  the  f  th  pixel  in  the  scan  line. 

n-io-  Ri 

n  -5  ~Ri 

X,  =  fi-Ri 
n+s~Ri 
h  +  io  ~  R, 

where: 

_  H-10  +  h-5  +  ?i  +  ?i  +  5  +  H  +  10 

R, - 5 

Since  /?,  is  a  local  average,  then  the  characteristic  increase  in  range  at  the  end  of  the  scan  lines 
has  a  decreased  effect.  The  feature  vector  X,  will  still  exhibit  some  slope  if  the  sensor  was  view¬ 
ing  a  roughly  horizontal  surface,  but  X,  will  be  centered  around  zero. 

There  will  be  some  holes  in  the  region  labeled  as  “target”,  pixels  where  the  Parzen 
classifier  chose  “non-target”.  Similarly,  there  will  be  isolated  points  in  the  background,  or  par¬ 
ticularly  likely,  in  the  clutter  regions,  mistakenly  labeled  as  “target”.  A  morphological  filter 
very  quickly  removes  these  single-point  (or  few-point)  errors,  yielding  a  binary  image  image 
labeled  as  “target”  in  regions  overlying  points  where  the  vehicles  were  originally  located,  and 
“non-target”  elsewhere.  Since  the  correct  “target”  regions  a^e  large  blobs,  and  the  false¬ 
positive  points  are  isolated,  it  is  a  simple  matter  to  simply  ignore  small  regions. 

Figure  4  shows  simulated  Topcat  data  for  a  cluttered  scene,  containing  three  T-72's  and  six 
trees.  Range  is  shown  as  intensity;  the  darker  splotches  are  the  tall  trees.  Careful  examination 
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of  the  target  regions  shows  the  turret  protruding  above  the  deck,  thus  slightly  darker,  and  the  bar¬ 
rel  projecting  out  the  front. 

Figure  5  shows  the  points  in  this  scene  initially  labeled  as  “target”  by  the  Parzen  estima¬ 
tor,  and  Figure  6  shows  the  detected  target  points  output  by  the  morphological  operator.  Notice 
how  the  “target”  regions  accurately  correspond  to  the  tanks  in  Figure  4. 

The  proposed  Topcat  system  would  be  in  a  fast-moving  missile  with  a  slight  look-ahead 
angle.  Figures  7  and  8  show  that  only  a  few  scan  lines  are  needed  to  detect  a  target.  In  this  way, 
Topcat  could  detect  a  target  and  still  have  adequate  time  to  fire  the  submunition  before  passing 
over  the  vehicle. 


CONCLUSIONS 

We  are  very  pleased  with  the  promising  results  we  have  shown.  This  algorithm  can  detect 
armored  vehicles  with  a  very  high  probability  of  detection  and  a  very  low  false  alarm  rate,  even 
when  the  vehicles  are  in  a  cluttered  scene.  Even  in  the  presense  of  clutter,  only  narrow  strips  of 
data  corresponding  to  a  few  scan  lines  are  needed  for  recognition,  meaning  that  the  algorithm 
would  be  amenable  to  real-time  application.  In  all,  we  feel  that  this  demonstrates  a  useful  Top 
cat  algorithm. 


Sequence  of  Ladar  Scan  Lines 


Figure  1 


Simulated  sensor  data  for  topcat  images 
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Figure  3 


The  feature  vectors  are  used  for  statistical  pattern 
classification.  The  Parzen  method  is  used  to  estimate  the 
probability  density  functions  for  target  and  non-target 
distributions. 


Overview  of  method: 


Comparison 
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Scene  0 


Simulated  TOPCAT  Image 


Missile  Speed:  274.3  m/sec 

Altitude:  22.5  m 

Spin  Rate:  150  rev/sec 

Sensor  Pulse  Rate:  70  KI  Iz 

Laser  Footprint:  6.4  mrad  x  1 .0  mrad 

Lookahead  Angle:  12° 

Field  of  View:  60° 
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Figure  4 


Figure  5 
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Detected  Target  Points 
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Simulated  TOPCAT  Image 
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Scene  2 


Simulated  TOPCAT  Image 
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